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We investigate the temperature dependence of the finite-field magnetization of the S = 1/2 bond- 
alternating XY model in the Random magnetic field along the z-direction having the Lorentzian 
distribution. The random-averaged free energy can be exactly calculated, which enables us to obtain 
the thermodynamical quantities. The temperature dependence of the finite-field magnetization shows 
various behaviors depending on the parameters 6 and T, where 6 and T are the bond-alternation 
parameter and the Lorentzian distribution parameter, respectively. 
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1. Introduction 

Many one-dimensional antiferromagnets have singlet ground states with the energy gaps in the 
excitation spectra. A typical example is the simple 5 = 1 chain having the Haldane ground state [1,2]. 
Another example is the bond- alternating 5 = 1/2 chain having the dimer ground state [3]. When the 
magnetic field applied to the spin-gapped system is so strong that the Zeeman energy exceeds the 
energy gap between the singlet ground state and the lowest state with finite magnetization, the system 
becomes gapless and can be described by the Tomonaga-Luttinger liquid theory. In this case, the low 
energy properties of the system are thought to be very similar to those of simple antiferromagnetic 
5 = 1/2 chain in zero magnetic field which is gapless. Thus, one may think that the temperature 
dependence of the magnetization may be very similar to that of the Bonner-Fisher's magnetic suscep- 
tibility [4]. However, as far as the temperature dependence of magnetization is concerned, the above 
naive expectation does not hold, as explained in the following. 

Hida, Imada and Ishikawa [5] studied the quantum sine-Gordon model with finite winding num- 
ber which can be related to the magnetization of spin chains. The 5 = 1/2 two-leg ladder was 
numerically studied by Wang and Yu [6] by use of the transfer-matrix renormalization group method, 
and also by Wessel, Olshanii and Haas [7] by use of the quantum Monte Carlo method. Both groups 
found a minimum-maximum (Min-Max) behavior (similar to that of Fig. 2(a)) of the magnetization as 
a function of the temperature, which is different from the Bonner-Fisher's behavior. Honda et al. [8] 
measured the magnetization of Ni(CsHi4N2)2N3(PF6) (often abbreviated as NDMAP) and observed 
the Min-Max behavior. Maeda, Hotta and Oshikawa (MHO) [9] performed the quantum Monte Carlo 
calculation for the 5 = 1 spin chain and found the Min-Max behavior. MHO also clarified the mech- 
anism and the universal nature of the Min-Max behavior. After MHO, the Min-Max behavior was 
found in several cases [10-13]. The Min-Max behavior was also found in the case that the origin of 
the spin gap is the spontaneous symmetry breaking [14]. 

Thus, the Min-Max behavior of the one-dimensional spin-gapped systems seems to be estab- 
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lished. However, the effect of the randomness on the Min-Max behavior has not been studied so far. 
Since the Min-Max behavior is based on the subtle balance between the distribution function and 
the density of states (DOS), as was explained by MHO, we should treat the randomness effects very 
carefully In this paper, we study the effect of randomness on the Min-Max behavior of the bond- 
alternating 5 = 1/2 XY model in the random magnetic field with the Lorentzian distribution which is 
exactly solvable. 



2. Model 



We investigate the model 

2N 2N 
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where 5 is the spin- 1/2 operator, / > is the coupling constant between neighboring spins, 27V is 
the number of spins, 8 is the bond-alternation parameter (0 < 8 < 1) and H is the uniform magnetic 
field along the z direction. The quantity Hj is the random magnetic field along the z direction with 
the Lorentzian distribution 

n H 2 + Y 2 

J 

We suppose that there is no correlation between Hi and Hj when i + j. Hereafter we set / = 1 (unit 
energy) and also - \,% - \ and gps - 1. 

The exact solvability of this model in the 8 - 0, Y = case was shown by Lieb, Schultz and 
Mattis [15], and that in the 8 > 0, Y = case by Pincus [16]. Nishimori [17] exactly solved the 
8-0, Y > case by use of the Lloyd method [18]. The full case 8 > 0, Y > was exactly solved by 
the present author [19]. 

3. Behavior of the Magnetization as a Function of Temperature 

Here we summarize the results of my previous paper [19]. The random-averaged free energy per 
one spin can be exactly calculated as 

f r°° i h 

F = — = -T dcop H (co) -r + -, (3) 
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where T is the temperature. The quantity p//(a>) is the DOS of Jordan-Wigner fermions given by 

Ph(oj) = pq(oj h ), a> H = (^ + H, (4) 

with 

i e + w 

Po(^) = — sin — — , (5) 

where A, B, 6 and (p are defined by 
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respectively. As is shown in Fig.l, the DOS po(ca) is 
everywhere non-zero when T > 0. In the absence of 
the random field (r = 0), we obtain 
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a> 



^_ — . = (tf<M<l) 

n V(^ 2 - 5 2 )(1 - CO 1 ) 
(otherwise) . 

(8) 

On the other hand, in the absence of the bond alter- 
nation (6 = 0), the DOS po(to) is reduced to Nishi- 
mori's result [17]. 

Once F is known, it is easy to calculate the mag- 
netization in the z-direction per one spin as 
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Fig. 1. DOS of the Jordan-Wigner fermions as a 
function of to for the r = 0, 0, 0.5, 1 .0 cases when 
6 = 0.5 and H = 0. 



Typical behaviors of M as functions of T are shown in Fig.2. Figure 3 shows the behavior patterns of 
M on the T - H plane which was obtained by the numerical analyses of eq.(9). 
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Fig. 2. Typical behaviors of the magnetization per one spin M = Mj 2N as functions of T for the (a) F = 0.0 1 , 
(b) F = 0.2 and (c) T = 1.0 cases when 5 = 0.5 and H = 0.6. The behavior of M in the Y = case is very 
similar to that in the F = 0.01 case. 



4. Discussion 



At low temperatures, by use of the Sommerfeld expansion [20], the magnetization M can be 
expanded as 
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where p' Q (H) - dpo(to) / dto^u and p["(H) = d i po{u))lda) i \ a>= H- Therefore the initial change of M is 
positive or negative according as po'(H) > or p~o'(H) < 0. As far as we take the terms up to T A in the 
rhs of eq.(10), the condition for the Min-Max behavior like Fig.2(a) is p' Q (H) < and p'q(H) > 0. 

When 5 = 0.5 and Y = 0, the Min-Max behavior of M is realized for 0.5 < H < 1 / V2. This behavior 
still holds even in the presence of small Y as shown in Fig. 3. Thus the shape of po(co) in Fig.l and the 
patterns of the behavior of M in Fig. 3 are consistent with each other. When T > 0, the excitation gap 
vanishes and, as a result, the divergence of the DOS also vanishes. Thus the existence of excitation 




gap and divergence of the DOS are not the necessary condition for the Min-Max behavior of M like 
Fig.2(a). Roughly speaking, the sharp double peaks of the DOS would rather be important. 

In conclusion, we have investigated the effect of the randomness on the Min-Max behavior for 
the first time. For the specific heat C we also observed interesting behaviors, for instance, the double 
peak behavior with a shoulder. Their details will be discussed elsewhere. 
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